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Abstract 

The rotation curve of the Large Magellanic Cloud, 
which we have derived from high-resolution HI 
position-velocity diagrams observed by Kim et al 
(1998), shows a steep central rise and flat rotation 
with a gradual rise toward the edge. Using the ro- 
tation curve, we have calculated the distribution 
of surface mass density, and show that the LMC 
has a dark compact bulge, an exponential disk, and 
a massive halo. The bulge is 1.2 kpc away from 
the center of the stellar bar, and is not associated 
with an optical counterpart. This indicates that 
the " dark bulge" has a large fraction of dark mat- 
ter, with an anomalously high mass-to-luminosity 
(M/L) ratio. On the contrary, the stellar bar has 
a smaller M/L ratio compared to the surrounding 
regions. 
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1 Introduction 

Rotation curves of dwarf galaxies are known to in- 
crease monotonically toward their outer edge, and 
the mass-to-luminosity ratio, and therefore, the 
dark matter fraction is larger than that in usual 
spiral galaxies (Persic et al 1996). However, little 
is known of their internal distribution of dynamical 
mass, because of the insufficient spatial resolution 
due to their small angular extents as well as for the 
low surface brightness in HI and Ha lines. 

HI observations of the Large Magellanic Cloud 
(LMC), the nearest dwarf galaxy, have provided us 



with a unique opportunity to study detailed kine- 
matics (e.g., Luks and Rohlfs 1992; Putman et al 
1998). Recently, Kim et al (1998) have obtained 
high-resolution HI kinematics from interferometer 
observations. They have shown that the rotation 
characteristics is nearly axisymmetric, and this fact 
allows us to derive a reliable rotation curve by using 
a position-velocity diagram along the major axis. 

In the present paper, we derive a new rotation 
curve for the LMC, and discuss the distribution of 
dynamical mass and dark matter. 

2 Position- Velocity Diagrams 
and Rotation Curve 

Fig. 1 shows the HI velocity field as reproduced 
from Kim et al (1998), superposed on a smoothed 
DSS (Digitized Sky Survey) image in B band. The 
cross indicates the kinematical center defined by 
Kim et al.. We have, here, determined a position, 
RA=5h 13.8m, Dec=-68°38' (J2000), as indicated 
by the asterisk, at which the velocity gradient at- 
tains the maximum, and around which the velocity 
field becomes most symmetric. From the symmetric 
velocity field, we may assume that the gas is circu- 
larly rotating around this kinematical center in the 
central 2 kpc radius region. The interferometer ob- 
servations may miss the intensities from extended 
gas disk. However, since the extended gas and the 
gaseous features detected with the interferometer 
are thought to be rotating at the same speed, we 
assume that the observed velocity field represents 
that of the total gas. 

In Fig. 2 we reproduce the PV diagram across 



a position close to the adopted kinematical center 
from Kim et al (1998). The thin line represents 
a rotation curve fitted by Kim et al. by using an 
AIPS utility, which adopts a simple function. How- 
ever, the innermost part of the PV diagram shows 
a much steeper variation from positive to negative 
velocities, which appears to be not well traced. In 
order to obtain a more precise rotation velocities in 
the innermost region, we apply the envelope-tracing 
method of maximum- velocities (Sofue 1996), which 
we have done by eye-estimates on the PV diagram. 
The newly determined rotation velocities are indi- 
cated by the thick curve, by which the steep veloc- 
ity variation is now well traced. Besides the bright 
features traced by this curve, a higher-velocity fea- 
ture is found at -60 to -70 km s _1 and +20'. If we 
trace this feature, the velocity variation around the 
center becomes steeper, and the rotation curve will 
have much steeper rise at the center. 

Using the thick curve in Fig. 2 and correcting for 
the inclination angle of 33°, we obtain a rotation 
curve as shown in Fig. 3. The outer rotation curves 
beyond the HI disk up to 8 kpc has been taken from 
Kunkel et al (1997). In the following discussion, we 
adopt the thick curve in Fig. 3, which is drawn 
by smoothing the observed curve (thin line). The 
thus-obtained rotation curve is found to be similar 
to those of usual disk galaxies (Rubin et al 1982; 
Sofue 1996, 1997; Sofue et al 1998, 1999) except for 
the absolute values. The rotation curve of the LMC 
is characterized by 

(1) a steep central rise, 

(2) a flat part in the disk, and 

(3) a gradual rise toward the edge. 

3 Surface Mass Distribution: 
Bulge, Disk, and Massive 
Halo 

Given a rotation curve from the center to the outer 
edge with a sufficient resolution, we can directly cal- 
culate the surface mass density (SMD)(Takamiya 
and Sofue 1999). This method is not intervened 
by any models such as the Plummer potentials and 
exponential disks. An extreme case is to assume 
a spherical symmetry: the rotation velocity is used 
to calculate the total mass involved within a radius, 
which is then used to calculate the SMD. Another 



extreme case is to assume a thin rotating disk: the 
SMD can be directly calculated by using the Pois- 
son equation (Binney and Tremaine 1987). We as- 
sume that the true mass distribution will lie in be- 
tween these two cases, which are indeed found to 
coincide within a factor of 1.5. Fig. 4 shows the 
calculated SMD for the disk (full line) and spher- 
ical cases (dashed line). The full line will better 
represent a disk component as the HI appearance 
indicates, while the spherical case would be better 
for an inner bulge. The obtained SMD distribution 
can be summarized as the following. 

(1) Dark bulge 

The distribution of SMD shows a compact, high- 
density component tightly concentrated around the 
kinematical center. This component is not associ- 
ated with an optical counterpart, and we call it a 
"dark bulge". The e- folding scale radius is about 
120 pc, when fitted with an exponential function. 
The mass involved within 240 pc radius, at which 
the SMD becomes equal to that of the disk com- 
ponent, is ~ 1.7 x 10 8 M Q . Except for its small 
radius, the dark bulge is as dense as a central bulge 
of normal disk galaxies: the SMD is estimated to be 
~ 1.0 x 10 3 -M Q pc at its scale radius, which is com- 
parable to or only by a factor of ten smaller than the 
values of usual galactic bulges (~ 10 3-4 M Q pc~ 2 ) at 
their scale radii (a few hundred pc) (Takamiya and 
Sofue 1999). 

(2) Exponential disk 

The SMD, then, is followed by an exponentially 
decreasing part with an e- folding scale radius of 1.0 
kpc. This "exponential disk" in mass is dominant 
up to a radius of 2 kpc, and extends to 3 to 4 kpc 
radius, being gradually replaced by an extended 
halo component. The total mass involved within 
2 kpc radius is ~ 2 x 10 9 M Q , consistent with the 
value obtained by Kim et al (1998). The SMD of ~ 
200M Q pc~ 2 at its scale radius is comparable to the 
values of usual disk galaxies (~ 300 — 5OOM pc~ 2 ) 
at their scale radii (^3 — 5 kpc). 

The total B magnitude of LMC is 0.6 mag. 
within an angular extent of 650' x 550', or within 
~ 4.3 kpc radius (de Vaucouleurs et al 1991). For a 
distance modulus of 18.50 (50 kpc) and B — V color 
of 0.55, we obtain a V magnitude —18.45. This 
leads to a V band luminosity of ~ 2 x 10 9 Lq. Since 
the mass involved within 4.3 kpc is ~ 4.5 x 10 9 M Q , 
we obtain a mean M/L ratio of about 2.2 for the 



disk component. 
(3) Massive halo 

The SMD in the outer halo obeys an exponen- 
tially decreasing function with an e-folding scale ra- 
dius of 5.2 kpc. The total mass involved within 8 
kpc radius is ~ 1 x 10 10 M Q . Since the outer enve- 
lope is not well visible in optical photographs (DSS), 
this may suggest again a high M/L ratio in the halo. 
This is consistent with the high M/L ratio in many 
other dwarf galaxies, for which gradually-rising ro- 
tation curves are generally observed (Persic et al 
1996). Kunkcl et al (1997) have suggested that the 
high velocities in the outer region could be due to a 
tidal effect with the Small Magellanic Cloud. Alter- 
natively tidal warping could also cause apparently 
higher (or lower) velocities, because the galaxy is 
nearly face on. In these cases, the presently esti- 
mated mass would be significantly overestimated. 

4 Discussion 

The dark bulge is significantly displaced from the 
stellar bar. No optical counterpart is visible, de- 
spite that the extinction is supposed to be not large 
for the nearly face-on orientation. Since the sur- 
face luminosity appears to be not significantly en- 
hanced from the average in the disk part, where 
M/L ratio is ~ 2, the bulge's M/L would be as 
large as - 20 - 5OM /L . Such a high M/L indi- 
cates a significant excess of the dark matter fraction 
in the bulge. As discussed in Section 2, if the faint 
and higher- velocity features in the PV diagrams are 
taken into account, the rotation curve will have a 
much steeper rise in the center than that traced by 
the thick line in Fig. 2. If this is the case, the 
bulge's mass density will be still higher, and so is 
the dark mass fraction. 

The center of the stellar bar is about 1.2 kpc away 
from the center of the dark bulge (Fig. 1). No sig- 
nificant mass enhancement, which might be related 
to the bar, is detected at radii 1 to 2 kpc in Fig. 
4. In fact, neither a stream motion in the velocity 
field, nor an anomaly in the PV diagrams are ob- 
served along the stellar bar (Kim et al 1998). We 
may argue that the M/L ratio in the stellar bar 
is significantly lower than that in the surrounding 
regions. 

We finally comment on an alternative possibil- 
ity to explain the displacement of the optical and 



kinematical centers. Ram-pressure stripping of the 
gas disk either due to an intergalactic wind (Sofuc 
1994) and/or a gaseous inflow from the Small Mag- 
ellanic Cloud (Gardiner et al 1994; Putman et al 
1998) could cause such a displacement. In order 
for the HI velocity field being kept unperturbed, 
it must occur as quickly as in a crossing time of 
the innermost disk, or within t ~ r/v rot ~2x 10 6 
years, where r ~ 120 pc is the radius of the in- 
nermost HI disk and i> ro t ~ 55 km s _1 is the rotat- 
ing velocity. This requires a wind velocity greater 
than i> ram ~ d/t ~ 500 km s , where d ~ 1 kpc 
is the displacement. Such a high velocity would 
be possible for an intergalactic wind. However, 
the ram-stripping condition, p ra m«ram > Pmv^ ot , 



seems to be not satisfied: for p r 
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cm - ", t> ram ~ 500 km s , phi ~ O.Itoh cm~ 3 , and 
v ro t ~ 55 km s _1 , we have 0.5 x 10~ 12 and 5 x 10~ 12 
dyne cm -2 for the first and second terms, respec- 
tively. Therefore, the ram-stripping appears not a 
likely origin. In any other circumstances, where no 
gravitational attraction is present, the inner rota- 
tion feature will disappear within about one million 
years. So, in so far as the velocity field is assumed 
to represent rotation, a high-density mass concen- 
tration at the kinematical center, and therefore, the 
presence of a dark bulge, will be an inevitable con- 
sequence. 

Another possible idea is to attribute the velocity 
gradients in the kinematical center to local velocity 
anomalies such as due to wiggles by HI shells and/or 
turbulence. However, Kim et al's modulus map of 
velocity residuals indicates that density waves and 
warping are more dominant to cause the residuals. 
Neither the largest HI shells nor 30 Dor are associ- 
ated with the largest velocity residuals and wiggles. 
Moreover, there appear no particular HI shells or 
associated starforming regions around the kinemat- 
ical center where the velocity variation is steepest. 

Hence, we here rely more on our assumption that 
the HI kinematics of the LMC directly manifests 
the dynamics of the galaxy. We also comment that, 
if we adopt the smooth rotation curve of Kim et 
al. (1998) as an alternative case, the mass distri- 
bution would be similar to the one represented by 
the disk component in Fig. 4. Even in this case the 
mass center is still significantly displaced from the 
optical bar, and no particular optical enhancement 
corresponding to the exponential disk is recognized 



around the kincmatical center. 

The author thanks T. Takamiya for calculating 
the SMD from the rotation curve data. He is also 
indebted to Prof. M. Fujimoto of Nagoya University 
for invaluable discussion. 
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position in this paper. Note the significant displace- 
ment of the bar from the kincmatical center. 

Fig. 2: Position- velocity diagram of the LMC along 
the major axis (Kim et al 1998). We trace the 
maximum- velocity envelope by the thick line to de- 
rive the rotation curve. 

Fig. 3: Rotation curve of the LMC. Note the steep 
central rise, and flat disk rotation. The thin line 
denotes the observation, and the thick line is a 
smoothed rotation curve, which wc adopt for cal- 
culating the surface-mass distribution. 

Fig. 4: Radial distributions of the surface mass den- 
sity for a thin disk assumption (full line) and for a 
spherical assumption (dashed line). There appear 
three components: (1) a dense and compact bulge 
at radii < 240 pc, (2) an exponential disk until 2 to 
3 kpc, and (3) a massive halo. 



Figure Captions 



Fig. 1: HI velocity field of LMC (Kim et al 1998) 
superposed on a smoothed DSS image in B band. 
The cross indicates the kincmatical center of Kim 
et al (1998), and the asterisk is the adopted center 
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